18.05 Lecture 14 
March 7, 2005 



Linear transformations of random vectors: Y — r{X) 



yi 




Xl 




= A 




Vn 




3^n 



A - n by n matrix, X = A'^ Y if det A ^ ^ A'^ = B 
Xl = biUi + ... + hinVn 

6ll... bin 



J = Jacobian = det 



bni-.. bn 



where b'^s are partial derivatives of Sj with respect to y. 



det B = det A'^ = -A- 
p.d.i. of Y: 



giy) 



\detA\ 



Example: X = {xi,X2) with p.d.f : 

f{xi,X2) = {cxiX2,0 < a;i < 1, < a;2 < 1; otherwise} 
To make integral equal 1, c = 4. 



Yi=Xi+2X2,Y2=2Xi+X2;A-. 
Calculate the inverse functions: 



1 2 

2 1 



det{A) = -3 



New joint function: 



^1 = -liYi - 2Y2),X2 = -\{Y2 - 2Fi) 



5(2/1,2/2) = X 4(-^(2/i - 2y2))(-^(t/2 - 2yi)) 
for < - 22/2) < 1 and < -^(1/2 - 2yi) < 1; 



0, otherwise} 



Simphfied: 



/(yi,2/2) = {^{yi - 2y2)(2/2 - 2yi) for - 3 < j/i - 2t/2 < 0, -3 < - < 0; 



0, otherwise} 
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X 



Linear transformation distorts the graph from a square to a paraUelogram. 



Note: From Lecture 13, when min() and max() functions were introduced, such functions 
describe engines in series (min) and parallel (max). 

When in series, the length of time a device will function is equal to the minimum life 
in all the engines (weakest link). 

When in parallel, this is avoided as a device can function as long as one engine functions. 

Review of Problems from PSet 4 for the upcoming exam: (see solutions for more details) 

Problem 1 - f{x) = {ce~^^ for x > 0;0 otherwise} 

Find c by integrating over the range and setting equal to 1: 



1= / ce-^'^dx 



Jo 




c = 2 



P(l < X < 2) = 2e-2^da; = e'^ - e"" 



Problem 3 - X - C/[0, 5],y = 0ifX<l;y = Xifl<X<3;y = 5if3<X<5 
Draw the c.d.f. of Y, showing ¥{Y < y) 



Y 



5 




x: 

Graph of Y vs. X, not the c.d.f. 



Write in terms of X ^ P(X-?) 



I 




1 3 5 
Cumulative Distribution Function 
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Cases: 



y < P(r <y) = P(0) = 

< 2/ < 1 ^ P(y < y) = P(0 < X < 1) = 1/5 

1< y < 3 ^ P(F < y) = P(X < y) = y/5 

3 < y < 5 ^ P(r < y) = P(X < 3) = 3/5 

y > 5 ^ P(y < 5) = ¥{X > 5) = 1 

These values over X from to oo give its c.d.f. 

Problem 8-0<a;<3,0<y<4 
c.d.f. F{x,y) = jiga;y(x2 +y) 

P(l < a; < 2, 1 < y < 2) = F{2, 2) - F{2, 1) - F(l, 2) + F(l, 1) 



or, you can find the p.d.f. and integrate (more complicated): 
c.d.f. of Y: P(y < y) = V{X < oo,Y < y) = V{X <'S,Y <y) 
(based on the domain of the joint c.d.f.) 

ny < y) = Tfe3y(9 + y) for < y < 4 

Must also mention: y < 0, V{Y < y) = 0; y > 4, P{Y < y) = 1 
Find the joint p.d.f. of x and y: 




J 



1 



1; 



Rectangle probability algorithm. 



f{x,y) = 



d^F{x,y) 
dxdy 



= {T^(3a;2 + 2y), 0<a;<3,0<y<4;0 otherwise} 




3 




93 
~ 208 



End of Lecture 14 
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